Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1 1348 1750 (N43W48 L= 82)   458    108        0 (1327/1712-1375/1773)  14.0 W43-W51 N41-N44
      2 1234 1402 (N25W22 L= 56)   351     86        2 (1198/1368-1270/1434)  15.9 W18-W25 N23-N25
      3  860 1377 (N31E03 L= 31)   297    103        0 ( 808/1362- 916/1389)  -9.5 E06-W00 N29-N32
      4  353 1349 (N38E45 L=350)   987    207        0 ( 294/1302- 423/1385) -17.1 E50-E37 N34-N41
      5  409 1094 (N21E40 L=355)   649    233        1 ( 315/1054- 490/1120) -19.7 E47-E33 N17-N24
      6  681  984 (N11E21 L= 13)   340     67        1 ( 644/ 966- 721/1007) -26.6 E23-E19 N09-N12
      7  223  917 (N13E59 L=335)   608    105        1 ( 199/ 875- 255/ 957) -68.2 E61-E55 N10-N15
      8  493  426 (S21E51 L=343)   433    120        1 ( 476/ 379- 519/ 472)  74.1 E52-E49 S24-S17
      9 1416  384 (S44W15 L= 49)   602    123        1 (1358/ 369-1470/ 411) -26.6 W09-W19 S47-S41
  

    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