Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1 1024 1511 (N36W02 L=121)   942    284        1 ( 920/1427-1189/1586) -33.7 E05-W13 N29-N41
      2 1458 1334 (N21W31 L=149)   450    110        0 (1413/1296-1509/1372)  26.6 W27-W35 N19-N23
      3  510 1336 (N25E35 L= 83)   728    144        2 ( 446/1312- 566/1367)  18.4 E40-E30 N23-N27
      4 1320 1311 (N21W21 L=139)   281     45        0 (1294/1286-1348/1334) -45.0 W19-W22 N19-N22
      5 1693 1274 (N16W49 L=167)   855    144        2 (1653/1222-1727/1321)  36.9 W44-W52 N12-N18
      6 1686 1111 (N06W46 L=164)   508     64        1 (1662/1073-1709/1144)  76.0 W43-W47 N03-N07
      7 1441  607 (S24W27 L=145)  1036    228        3 (1388/ 476-1495/ 732)  73.6 W24-W30 S33-S16
      8  543  621 (S20E35 L= 83)   256     66        0 ( 510/ 592- 561/ 648) -63.4 E36-E33 S21-S17
      9  744  406 (S36E25 L= 93)   918    270        2 ( 671/ 356- 824/ 477)  28.1 E32-E17 S39-S30
  

    OPENR, u, binfile, /COMP, /GET_LUN
    head = FLTARR(6) ;[xsize,ysize,radius,hh,mm,ss]
    READU, u, head
    data = INTARR(FIX(head[0]),FIX(head[1]))
    READU, u, data
    FREE_LUN, u
    

• each filament has the value of Fil_ID
• all values are in pixel

  c = WHERE(data EQ Fil_ID, area)
  f = '(2(A,I6))'
  PRINT, FORM=f, "Area of Fil_ID", Fil_ID," =", area 

  data = ROTATE(data,7)
  Solar_y = 2048-Solar_y