Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1 1418 1567 (N32W39 L= 33)   544    207        0 (1338/1547-1490/1600)  -4.4 W33-W46 N30-N33
      2 1113 1443 (N30W13 L=  7)   757    152        2 (1050/1399-1195/1479)   9.5 W08-W20 N28-N32
      3 1676 1359 (N14W51 L= 44)  1318    191        3 (1605/1296-1750/1409) -48.0 W45-W56 N08-N18
      4  536 1075 (N16E31 L=323)   522    144        0 ( 473/1046- 597/1115) -24.0 E35-E26 N13-N19
      5 1162  909 (S03W06 L=360)  1795    418        7 (1113/ 727-1208/1104)  79.7 W02-W09 S13-N09
      6 1457  908 (S08W25 L= 18)   707    138        3 (1393/ 850-1541/ 957) -48.4 W20-W30 S12-S03
      7 1326  796 (S12W15 L=  8)   446    101        0 (1303/ 742-1353/ 854)  66.8 W12-W17 S15-S08
      8 1336  594 (S26W13 L=  7)   575    104        0 (1308/ 544-1367/ 653)  50.2 W11-W16 S28-S22
      9  983  351 (S37E16 L=337)  1225    301        2 ( 908/ 279-1031/ 468)  43.0 E24-E09 S42-S28
  

    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