Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1 1403 1378 (N24W28 L=342)   648    139        2 (1328/1342-1469/1409) -24.0 W23-W32 N21-N26
      2 1581 1299 (N18W40 L=354)  1182    224        4 (1517/1235-1639/1370) -52.1 W35-W43 N13-N22
      3 1835 1278 (N14W65 L= 18)   969    105        3 (1802/1236-1867/1328) -63.4 W60-W68 N10-N16
      4  266 1108 (N11E54 L=259)  1134    220        4 ( 217/1068- 311/1179)  66.8 E59-E49 N08-N15
      5  879 1122 (N11E09 L=305)   219     63        0 ( 853/1101- 903/1155)  45.0 E10-E06 N10-N13
      6  531 1054 (N08E32 L=282)   259     39        0 ( 506/1037- 549/1075)  33.7 E33-E30 N07-N09
      7 1865  937 (S08W64 L= 18)  1228    210        2 (1810/ 893-1911/ 970) -23.2 W56-W71 S10-S04
      8 1520  753 (S16W32 L=346)  1112    293        1 (1410/ 639-1619/ 892)  47.0 W24-W38 S22-S07
      9  161  819 (S07E71 L=243)   614    140        1 ( 140/ 802- 188/ 845) -15.9 E74-E66 S08-S05
     10 1610  554 (S30W42 L=355)  1364    212        2 (1511/ 513-1684/ 597)  12.1 W33-W48 S32-S27
     11  424  542 (S25E47 L=267)   748    145        2 ( 371/ 502- 477/ 581) -26.6 E51-E43 S27-S21
  

    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